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Preface 

In  modern  planning,  predicted  performance  of  programs,  integrated  operations, 
and  individual  pieces  of  equipment  play  an  important  role.  This  is  particularly  true 
of  those  defense  items  associated  with  the  general  problem  of  hitting  a  target. 

The  table  of  salvo  kill  probabilities  for  square  targets  presented  herein  was  pre¬ 
pared  by  the  National  Bureau  of  Standards  Institute  for  Numerical  Analysis,  Los 
Angeles,  at  the  request  of  North  American  Aviation,  Inc.,  and  with  the  support  of  the 
Flight  Research  Laboratory,  Air  Research  and  Development  Command,  United  States 
Air  Force.  However,  U  should  be  useful  to  other  organizations  engaged  in  such 
operations  as  operation  analysis,  strategic  planning,  engineering,  and  field-test  activities, 
as  well  as  missile-development  programs.  Hence,  following  the  initial  preparation 
of  a  small  number  of  copies  to  meet  the  needs  of  the  immediate  requester,  formal 
publication  seemed  warranted  to  make  the  results  more  widely  available.  Although 
the  terminology  and  notations  used  are  specifically  directed  to  the  defense  problem, 
the  results  are  not  limited  to  that  field. 

A.  V.  Asrm,  Director. 

May  28,  1Q54.| 
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1.  Iiitrodiietioii 

].l  IMathcniatical  Formulation 


In  the  prnhiein  of  niiniii^  n  \von|)on  to  hit  n  (ni‘}'et,  two  oiroi*s,  nssunieil  to  he  iiidepetuleiit,  nro 
consiiloml.  One  is  tiio  dispersion  of  the  missile  firetl.  It  is  a  function  of  the  hnllistic  charnct eristics  of 
the  particular  missile  and  of  the  nieteoroloi;ieal  eomlitions  existing  at  the  time  and  place  of  release.  It 
may  well  be  considered  ns  a  random  error.  In  this  treatment  of  the  subject,  the  x  and  y  errors  nro 
assumed  to  bo  independently’  normally  distributed  with  center  at  the  point  of  aim  and  a  standard 
deviation  of  and  in  /  and  y  directions,  respectively. 

The  aiming  error  is  the  second  error  affecting  the  residts.  The  x  ami  y  com|M)nonts  of  the  aiming 
error  also  are  assumed  to  he  indcpeiulently  normally  distributed  with  renter  at  (x„,y^  and  standard 
deviation  of  ^,4^  and  ^4^,  respectively.  For  a  particular  salvo,  the  effect  of  this  aiming  error  is  definitely 
biased  in  some  direction.  Hence  we  may  refer  to  this  error  as  a  systematic  error. 

Assume  a  sipiare  taigct  of  side  2a  whose  center  is  at  the  origin  and  whose  siiles  are  parallel  to  the 
coordinate  u.\es.  In  this  case  the  probability  of  hitting  the  target  by  a  missile  aimed  at  is 


(1) 


Ia'I  I’k  be  the  probability  that  one  missile  destroys  the  target  or  makes  it  useless.  If  N  missiles 
aimed  at  are  fired  in  salvo  and  if  they  act  independently,  the  probability  that  the  target  will  be 
destroye*!  is 


This  will  be  calleil  the  salvo  kill  probability  of  missiles  ainie<l  at  ({,77). 

In  practice,  one  does  not  know  the  true  aim  point  ({,77).  Wc  have  assumed  the  aiming  distribution 
to  be  normal  with  center  at  (ro,7/o).  Therefore,  the  salvo  kill  probability  of  A’ missiles  released  in  salvo 
is,  under  the  above  assumptions, 


PsK-= 


1 _ 


((-».)■  1  (n- 


J 


(2) 


The  (quantity  Psk  «  function  of  the  paramctei’s  Xq,  ija,  VAy,  ^a,,  and  A’. 

In  practical  applications  one  seldom  has  targets  of  fi.\cd  shapes,  i.  e.,  squares,  rectangles,  circles, 
oblongs,  etc.  Hence  the  choice  of  target  shape  is  quite  arbitrary.  The  orientation  of  the  target  relative 
to  the  center  of  the  aiming  point  will  certainly  affect  the  result.  However,  it  is  believed  that  conclusions 
based  on  application  of  these  probabilities  will  bo  insensitive  to  the  precise  orientation  of  the  target 
relative  to  the  point  of  aim.  Consequently,  we  have  assumed  the  square  target  is  always  oriented  so 
that  one  of  its  sides  is  normal  to  the  line  joining  the  center  of  the  target  with  the  center  of  the  aiming 
distribution.  Tims  Xo  is  taken  to  be  zero  in  the  tabidation  presented  herein,  also,  the  disi''ibutions  were 
assumed  to  have  equal  variances  for  x  and  y,  i.  e., 


VA=-s«Ay-^«Ay=9Ayy2=OAyy^  y 
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1.2  Numerical  Analysis 

Under  the  bIjovo  nssiimptionfi  I  lie  expremion  for  /’«((,  n)  «**  given  in  (1)  can  be  simplified.  Define 


Then  /'*({, i;)  may  he  written  ns 


Aj)=  I  < 

>  2wJo 

itten  as 


Numerical  integration  of  /V  (e(|  H)  111H3’  be  accomplished  as  follows: 

A  grid  of  squares,  each  side  of  which  is  a/n,  is  superimposed  on  the  coordinate  system  as  show'n  in 
figure  1.  Denote  by  the  center  of  square  whose  vertices  are  {{iainjain),  ((1  + !)<*/«, ja/n), 

((.i+\)aln,(j+l)aln),  and  ({aln,(j+l)aln). 

Ijct  w  be  chosen  such  that  (J(tri)  may  he  approximated  by  0((itVj)  for  ta/n  <(<({+ l)a/n  and 
ja/n  <  i)<  0'+ !)(//».  This  assumes  that  the  grid  is  sufliciently  small,  for  the  accuracy  aimed  at,  to  con¬ 
sider  the  salvo  kill  prohability  constant  within  the  square  Ttj.  With  this  assumption  the  integral  in 
(3)  may  be  approximated  by  a  double  sum 

The  term  in  the  hiockots  represents  the  probability  that  the  point  of  aim  will  lie  in  grid  T,j. 

This  term  denoted  by  j)  may  bo  written  as 

Toking  note  of  the  existing  s.vnunetry,  we  obtain 


Q(i,ri,)l\(hjr  (9) 


-1  -f  -I  .  Its 
1*0 
(a) 

Fiocaa  1. 
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In  the  associated  tables  the  following  conditions  were  assumed:  (1)  Tlie  random  and  systematic  errors 
are  each  circular  normal  distributions  with  equal  variances,  and  (2)  the  center  of  the  aiming  distribution 
is  on  the  (i.  e.,  at  (0,ya)>  Making  use  of  these  assum|itioiis,  the  c|unntities  required  are 


^«(fVciyoie,4ie«iA)*2  2  £ 


These  equations  a’ere  evulunte<l  for  the  follow  ing  vnhuw  of  the  pnniineters: 

/V=0.1,  0.4,  0.7,  I.O 
y„=0,  a,  2rt,  4rt,  7a,  llo,  KWi,  22a 
9n=a,  2a,  4a,  7a,  I  to,  IG<i,  22a 
€A—a,  2a,  4o,  7a,  I  to,  16a,  22a 
A’=  1,  5,  10,  25,  50,  100,  150,  200. 


(10) 


1.3.  Method  of  Computation 

In  order  to  solve  cq  (10)  by  modern  high-speetl  computers,  such  ns  SWAC  (Standards  Western 
Automatic  Computer),  certain  decisions  and  'or  npproximntions  must  be  made.  The  quantities  PM((t,Vj) 
and  /\(<,  j)  may  be  obtained  from  tabulated  values  of /(r).  Such  n  method  is  inefTicient  for  SWAC 
solution.  Consequently,  a  Clicbyshev  |)olynominl  e.xpnnsion  was  usctl  to  nppro.\imate  /(»).  This 
polynomial  is 

/'(/)=Oo+a,H«a  {’+.  .  .+0.,?'*;  {=(/-3)/3 
and 

a„= 0.4980  5021  flT=-.2131  8912 

a,  =  0.0 133  0987  «»=0.3101  6832 

aj= -0.0598  4802  «,=  - 0.0541  0048 

03=0.1590  0224  fl,«=-0.1812  5824 

04= -0.2087  2400  fl„=0.0980  5824 

Os= 0.2751  8080  o„= 0.0309  0490 

o«  =-0.0858  9280  o„=- 0.0282  0240 

This  function  differs  from/(j-)  by  a  small  amount.  The  error /(/)—/'(/)  is  shown  in  figure  2. 


♦  »o  I  to*» 

•  ♦  I  !•*» 

)  0 
m 

•  -  B  !•*» 


Fiovaa  2.  Graph  of  /(x)-/'(*). 
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TliP  grill  sizp  phown  so  ihnf  for  the  npfiinu*y  <li*Hirp«l,  rotihl  Ik*  cniisidorcd  cnnslnnt  williiu 

(•ncli  mrnil>rr  of  the  grid.  TIliA  wns  accotnplislipcl  ii)*  Holrrting  n  grid  siirh  that  tlip  vnliip  of  Psk  ^vas  not 
pimngpd  (to  tliP  nccurnpy  dpsirpd)  when  tho  grid  8i*p  was  doiddpd.  Analysis  of  this  inptho«l  loails  to  tho 
poni'lusion  that  tliP  grid  sizp  depends  primarily  on  Xumerienl  ooinputation  shows  that  the  following 
ehoiee  of  n  woidd  yield  five  decimal  results,  which  would  l)e  unchanged  by  doidding  n. 


a 

2a 

4a 

7a 

11a 

IGa 

22a  1 

•  1 

n 

8 

4 

3 

2 

1 

B 

'A 

1 

Truneation  is  always  a  problem  in  numerical  compulation.  In  this  problem  the  maximum  con¬ 
tribution  to  the  integral  naturally  comes  from  the  center  of  the  combined  distribution.  Tliis  center  is 
located  l>etween  the  origin  and  (O.^o)  and  depends  on  N,  the  number  of  missiles  in  the  salvo.  In  tho 
machine  solution,  SWAC  was  coded  to  search  this  area  and  the  surrounding  area  until  the  contribution 
from  each  member  of  the  grid  fell  below  a  specified  limit.  Tliis  limit  was  determined  by  computing  Pmk 
for  cases  known  to  converge  slowlv,  using  a  trial  limit.  This  trial  limit  was  halved  until  the  difference  of 
two  successive  values  of  Psk  was  negligible.  In  this  manner  the  limit  of  10~*  was  determined. 

Cliecking  techniques  are  a  must  in  numerical  computation.  Elementary  analysis  applied  to  tho 
expression  for  Pax  (eq.  3)  reduces  the  double  integral  to 


j^ji^(-ro+a)V2)  q)%'2^J|y^(ife+a)V2^ 


where  +  and  A'=l.  Thus,  the  case  A’^=l  can  be  obtained  directly  from  tubulated  values  of 
/(/).  This  provides  an  excellent  independent  chwk  for  the  SWAC  solution. 

Further  analysis  shows  that  Pax  approaches  unity  ns  N  becomes  infinite.  This  gives  a  qualitative 
chock  on  the  results.  Another  cheek  was  applied.  Pax  f<*r  A^=l,  5,  10,  25,  50,  100,  150,  200  were 
obtained  simultaneously  by  forming  powers  of  [1  —/**/**)•  To  insure  that  the  routine  was  forming  tho 
correi^t  powers,  several  values  of  Pax  were  computed  by  hand  and  compared  with  the  corresponding 
ninchine  solution.  The  agreement  was  o.xceileiit. 

As  a  further  check,  the  salvo  kill  probabilities,  Pax,  were  eomputed  twice  by  SWAC.  The  agreement 
of  the  two  sets  of  results  represents  a  further  assurance  of  eorrect  computation.  As  a  final  check  cards 
were  punched  from  the  galley  proofs  and  compared  with  the  original  punched-card  tabic. 

Considering  the  methods  applied,  the  selection  of  the  grid,  the  truncation  criteria,  and  the  checks 
applied,  we  believe  the  tabular  entries  given  herewith  are  correct  to  three  decimals.  The  fourth  decimal 
is  in  error  by  at  most  two  units. 


1.4.  Choice  of  Parameters 

Many  users  will  perhaps  question  the  particular  choice  or  parameters.  A  word  of  explanation  is  in 
order.  In  tho  use  of  high-speed  computing  equipment,  especially  on  SWAC,  memory  storage  space  is 
at  n  premium.  In  order  to  have  as  much  storage  as  possible  the  parameters  were  generated  witliin 
SWAC  rather  than  stored.  This  lead  to  a  choice  of  parameters  that  can  lie  generated  by  elementary 
means. 

1.5.  Interpolation 

This  table  was  designed  to  meet  a  specific  purpose,  and  it  was  not  intended  to  provide  an  interpolable 
table  in  the  directions  of  ya,  vx,  or  xa.  Tlie  table  would  have  to  be  enlargetl  many-fold  to  make  it  inter¬ 
polable  in  all  these  directions.  However,  in  certain  regions,  it  is  possible  to  interpolate  in  the  Px  direc¬ 
tion,  for  fixed  values  of  the  other  parameters.  The  functions  for  N—  1  are,  of  course,  linear  in  Px.  For 
AT  greater  than  one,  the  functions  cease  to  be  linear  in  Px.  The  departure  from  linearity  increases  with  AT. 
Consequently,  the  accuraey  of  linear  intn'polation  will  decrease  with  increasing  N. 
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The  coniputntion  of  a  table  of  this  ningnitiide  represents  the  rooperntion  of  several  individuals, 
each  of  M’hom  deserves  niueh  credit.  E.  C.  Yowell  gave  consitlerahie  assistance  in  every  phase  of  the 
task,  Dan  Teichroea’  was  consulted  on  statistical  matters  pertinent  to  the  preparation  of  tho  table. 
The  coding  and  macliine  operation,  under  the  general  giddanco  of  Roselyn  Lipkis,  were  performed  by 
Nan  Reynolds  and  Shirley  Marks.  Hand  computation  of  the  numerous  check  values  was  performed  by 
Marvin  Howard  and  his  staff.  F.  Hollander  and  Patricia  Bremer  and  their  punched  card  machine  staff 
assisted  in  the  preparation  of  the  cards  and  the  listing  of  results.  Gertnule  Blanch  e.\amined  the  table 
for  interpolability  and  recommended  tho  appropriate  fonnat. 

Arn’old  D.  Hestenes. 

Los  Angeles,  Calif.,  April  2,  1053. 
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